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Strong Rashba spin-orbit coupling (SOC) of two-dimensional electron gas (2DEG) at the oxide 
interface LaAlOa/SrTiOa underlies a variety of exotic physics, but its nature is still under debate. 
We derive an effective Hamiltonian for the 2DEG at the oxide interface LaAlOa/SrTiOa and find a 
different anisotropic Rashba SOC for the d^z and dyz orbitals. This anisotropic Rashba SOC leads 
to anisotropic static spin susceptibilities and also distinctive behavior of the spin Hall conductivity. 
These unique spin responses may be used to determine the nature of the Rashba SOC experimentally 
and shed light on the orbital origin of the 2DEG. 


PACS numbers: 73.20.-r, 71.70.Ej, 72.25.Mk 

The discovery of a high mobility two-dimensional elec¬ 
tron gas (2DEG) at the interface between two band insu¬ 
lators LaAlOa and SrTiOa (LAO/STO) [1] has attracted 
increasing attention [2]. However, the origin of the 2DEG 
is still under active debate. According to the intrinsic po¬ 
lar catastrophe mechanism, there should be a half elec¬ 
tron (per unit cell) transfer from the top surface layer 
of LAO to the LAO/STO interface. The resulting car¬ 
rier density at the interface is roughly 3.5 x 10^“^ cm“^, 
which mainly comes from the three t 2 g orbitals of Ti in 
STO. Several transport experiments, however, estimate 
that the carrier density is only 10% of that due to the 
polar catastrophe mechanism [SHS]. In addition, it has 
been proposed that electrons in the dxy orbitals, which 
are confined in the xy plane, are more likely to become lo¬ 
calized at the interface due to the impurities or electron- 
phonon coupling, while those in the dxz and dyz orbitals 
are itinerant and contribute to transport [S] . Within this 
scenario, the localized and itinerant electrons would ac¬ 
count for the observed magnetic order [7] and supercon¬ 
ductivity EHII], respectively. It is therefore important to 
understand the transport properties of the 2DEG. Other 
mechanisms, such as oxygen vacancies and polar 

distortion [Ml [15], have also been proposed. 

Recent magnetotransport experiments have provided 
us insight into the 2DEG at the oxide interface. In par¬ 
ticular, a strong and field-tunable Rashba spin-orbit cou¬ 
pling (SOG) was observed [THIITTI and was modeled using 
the standard fc-linear form |18j . i.e.. 

Hr = Ak(A; x cr) ■ z . (1) 

Based on this fc-linear Rashba SOC, theoretical works 
have predicted a variety of unusual effects, such as Fulde- 
Ferrell-Larkin-Ovchinikov-type superconductivity coex¬ 
isting with ferromagnetism m, spiral magnetic order 
and skyrmions , and the spin Hall effect [23]. How¬ 

ever, a very recent magneto-conductivity measurement 
has suggested the possibility of a A:-cubed Rashba SOC 
of the 2DEG at the oxide interface |25l |26| . Accordingly, 
some authors proposed a fc-linear Rashba SOC for the 


dxy orbital EZllSH] and a /c-cubed one for the dxz and 
dyz orbitals [2H|. On the other hand, first-principles cal¬ 
culations combined with the envelope function method 
have found an anisotropic nonparabolic spin-split sub¬ 
band structure for the dxz and dyz orbitals [55], which 
could not be explained by the standard fc-cubed Rashba 
SOC. Thus, a detailed investigation of the low energy ef¬ 
fective model and the nature of the Rashba SOC is highly 
desirable. 

In this Rapid Communication, we present a detailed 
derivation of the effective Hamiltonian of the 2DEG 
at the oxide interface. We find a different anisotropic 
Rashba SOC of the following form, 

Ht « {kl - kl) (fc X <t) • z, (2) 

for the dxz and dyz orbitals, and a standard fc-linear 
Rashba SOC for the dxy orbital. The anisotropy of the 
Rashba SOC naturally leads to anisotropic spin suscepti¬ 
bilities that have been observed experimentally [Mill]. 
We also show that this anisotropic Rashba SOC results 
in different behavior of the spin Hall conductivity (SHC) 
when compared to the standard fc-linear and fc-cubed 
Rashba SOCs. These distinctive spin responses can be 
used for determining the nature of the Rashba SOC in 
experiments and to shed light on the orbital origin of the 
2DEG at the LAO/STO interface. 

We begin by constructing the low-energy effective 
model of the 2DEG at the LAO/STO interface around 
the r point in the Brillouin zone. The 2DEG is formed 
from the d orbitals of the transition-metal Ti. Here we 
focus on the three t 2 g orbitals, namely, dxy, dxz, and dyz, 
since the Cg orbitals are pushed up about 2 eV higher 
than the Og orbitals by the octahedral crystal field. On 
the xy plane, electrons in the dxy orbital can hop along 
either the x or y direction to the dxy orbitals on the neigh¬ 
boring Ti, while electrons in the dxz (dyz) orbital can hop 
to its neighbor only along the x (y) direction. Thus, the 
corresponding hopping Hamiltonian can be expressed in 
the following matrix form, 
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/ h{k) 0 0 

i?o = I 0 —2tcoskx 0 

y 0 0 —2t cosky 


( 3 ) 


where h{k) = —Ae — 2t{coskx + cosky), t = tp^/Apd is 
the effective hopping parameter between nearest neigh¬ 
boring Ti, Apd is the splitting between the oxygen p and 
Ti t 2 g energy levels, and is the difference in the on¬ 
site energies between the dxy orbital and the dy^/dxz or¬ 
bital. Note that since dxy is even, and dxz and dyz are 
odd under the operation z — >■ —z, hopping between these 
two sets of orbitals is prohibited in the presence of the 
mirror symmetry. 

To model the effect of the SOC, we introduce 
the atomic SOC • cr in the basis 
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where ^ denotes the strength of the atomic SOC. cr refers 
to the spin degree of freedom, while I is the orbital an¬ 
gular momentum of the electron. 

Finally, there is a mirror symmetry breaking at the 
interface due to the polar displacement of Sr and Ti 
atoms relative to the oxygen octahedra, which leads to 
the Rashba SOC. Physically, the mirror symmetry break¬ 
ing can induce the hopping process from the dxz (dyz) 
orbital to the dxy orbital via the Px (py) orbital of oxy¬ 
gen. The corresponding Hamiltonian can be written as 

[IZlllH] 


H^ = l 


0 — 2 isinfcy —2isuikx 

2i sin ky 0 0 

2i sin kx 0 0 


O cTo, (5) 


where 7 refers to the effective hopping amplitude between 
the dxy orbital and the dxz and dyz orbitals. (Tq is the 
2 x 2 unit matrix in the real spin space. 

The total tight-binding (TB) Hamiltonian including all 
three parts is given by 


Htb = FZo + • 

There are three pairs of degenerate bands at the T point, 
which are plotted in Fig. 1(a) using the parameters given 
in Ref. [2H1 It can be seen that the energy contour of the 
middle two bands has a strong anisotropy as shown in 
Fig. 1(b), whereas the lowest two bands are isotropic as 
shown in Fig. 1(c). Note that the splitting of the two 
lowest-energy bands due to the Rashba SOC is unnotice- 
able for the given energy. 


To derive the effective Hamiltonian, we apply the 
quasidegenerate perturbation theory [30]. Up to lead¬ 
ing order in the SOC strength we obtain the effective 
Hamiltonian for the top pair of bands. 


ddtop (^) — z Otop {k X cr) • z, (6) 

■^^top 

the middle pair of bands. 


^^mid (fe) = 




+ Omid {kl - kl) {kx cr) ■ z, (7) 


and the bottom pair of bands, 


Hbot (k) = - -abot {k X cr) ■ z, 

znihot 


( 8 ) 


where (mtop, mmid, Wbot) and (atop, ciniid, abot) are the 
effective masses and Rashba SOC strengths for the top, 
middle, and bottom pairs of bands, respectively (all their 
specific expressions are given in the Supplemental Ma¬ 
terial m)- The top pair of bands is a mixture of all 
three t 2 g orbitals. The bottom pair mainly comes from 
the dxy orbital. The middle pair is a hybridization be¬ 
tween the dxz orbital and the dyz orbital. It is also clear 
that the bottom pair of bands has the Zc-linear Rashba 
SOC that was proposed by previous works [171I2HI- This 
concentric isotropic Fermi contour of the dxy orbital had 
also been demonstrated at the surface of bare SrTiOa 
[32]. In the middle pair of bands, the Rashba SOC 
becomes anisotropic and has a Zc-cubed energy disper¬ 
sion |33j . Two recent angle-resolved photoemission ex¬ 
periments had already observed the anisotropic Fermi 
contour of the dxz orbital and the dyz orbital at a high 
carrier density [34l[35]. Note that the effective Hamilto¬ 
nian of each pair of bands is constructed with respect to 
its own bottom edge. 

The anisotropic Rashba SOC for the dxz and dyz or¬ 
bitals in Eq. Q is our main result. In the rest of this 
Rapid Communication, we will study its effects on the 
static spin susceptibility and the spin Hall conductiv¬ 
ity [5S1 [37]. For convenience, we redefine the corre¬ 
sponding effective mass m = rrimid/k?' and Rashba SOC 
strength /3 = ■ The effective Hamiltonian for the 

middle pair can be recast into 


^ {kl - kl) k_ 


Hniid{k)= I . 0 ^ 3 /, 2 "* 1,2 ^ 


(9) 

with k± = kx A iky. Some simple algebra leads to the 
eigenvalues of i7mid {k), 


^ks 


h^k^ 

2m 


+ sjdfi^k^ |cos 20 fc|, 


( 10 ) 


and the corresponding eigenvectors. 


4>ks = 


( 11 ) 
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FIG. 1. (Color online) (a) Band structure of TB model de¬ 
scribing the oxide interface. Energy contours near the F point 
for energies (b) Ejt = —2 and (c) E/t = —3. Parameters 
are adopted from Ref. 1281 = —0.56, = 0.035, 

7/t = 0.072. 


where the spinor is given by rjks = (—l)^ /a/S, 
s = ±1 is the chirality index, and is the area of 
the 2DEG with = cos20fc/|cos20fc| = ±1, 9^ = 
arctan(fcy/fca;). Since our model is only valid around the 
r point, we would like to introduce a momentum cutoff 
kc = l/imfiP via the turning point of the energy disper¬ 
sion Ek = k^fi^/2m — The corresponding energy 

of this turning point is given by etum = l/54m^/3^. 

In general, the free spin susceptibilities can be written 
as 


XijiQ) = cTjG {k + q,ujn)], 

n,k 

( 12 ) 

where ai are the Pauli matrices with i = x,y, z, G (fe, u!„) 
is the Matsubara Green’s function of an electron with 
momentum k and frequency and /r_B is the Bohr 
magneton. After carrying out the standard analytic con¬ 
tinuation and frequency summation (more details of the 
derivation can be found in the Appendix of Ref. I39p . we 
can find the static spin susceptibilities in the limit q —> 0, 
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FIG. 2. (Color online) The spin susceptibility of 2DEG with 
the anisotropic Rashba SOC as a function of the Fermi energy 
Ef in units of etum (measured from the bottom of the middle 
pair of bands). We set the dimensionless effective mass of the 
electron m = 1, the dimensionless Rashba SOC parameter 
P = 0.01, and the temperature T=5 K. 


All of the other components vanish due to the symme¬ 
try of Fermi surface. The out-of-plane component x^z is 
the so-called van Vleck susceptibility and originates from 
the virtual inter band transition. The in-plane compo¬ 
nent Xxx or Xyy contains both the intraband contribution 
(the Pauli susceptibility) and the interband contribution 
(the van Vleck susceptibility). 

Numerical calculations of the spin susceptibilities of 
2DEGs with the anisotropic Rashba SOC show two main 
features, as shown in Fig. 2. First, the spin suscepti¬ 
bilities are anisotropic, i.e., x^z 7 ^ Xxx- Secondly, the 
spin susceptibilities have a strong Fermi energy depen¬ 
dence. Note that the momentum cutoff kc is used in our 
numerical calculations. 

Previously the anisotropic spin susceptibility was also 
found using the fc-linear Rashba model [55] • However, 
the spin susceptibility is anisotropic when only the lower 
Rashba spin-split band is occupied. As soon as both spin- 
split bands are occupied, the spin susceptibility becomes 
isotropic |39j . As such, the anisotropy only shows up in a 
small energy window. In contrast, the spin susceptibility 
in our model is always anisotropic (up to the turning 
point when the model is no longer valid) |40|. Therefore, 
our result may provide an alternative explanation for the 
observed magnetic anisotropy nniiii]. 

Let us now turn to calculate the SHC of the 2DEG 
with the anisotropic Rashba SOG. The general spin con¬ 
ductivity tensor in the spin space is given as 


27rL2 




jCTi 


(16) 


where (k) = Ex (fc) — Ep is the energy of the electron 
measured relative to the Fermi energy Ep, and the su¬ 
perscript 0 indicates the spin susceptibility with q — 0. 


where is the vertex function in the spin 

space and Kx = G^JxG^ is the vertex function in the 
eigenvectors space of H,nid{k). G^ and G^ are the re- 
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tarded and advanced Green’s function of 2DEG, 


Born approximation can be written as 


qA 


(^) = 7 Z/ (e)= / 

^ ^ks ^ ^ks 


■ IT] 


(17) 


where 77 is a positive infinitesimal. Jq, stands for the 
velocity operator in the eigenvectors space and is given 
by Jcc = U^jaU, where ja = eva is the current operator 
of electron in the spin space and Va = dH^i^/d (hka) 
refers to the velocity operator with a = x,y. The spin 
current operators are represented by 

Ja , (18) 

where {A,B} = AB + BA is an anti-commutator, and 
the 2 x 2 unitary transformation matrix is of the form 


J_ / -zcfce KfcC \ 

1 1 ; 


(19) 


After taking the trace over the spin degree of freedom, 
we have the nonzero component of the intrinsic SHC as 


ebXhJ 

IGtt^ 


k^dk 


n2Tr 


Ep — h?k^/2m Jq 


Cfc sin 


X cos2ek[6{EF - Ek-) - S{Ef - ek+)]d9k, (20) 


which indicates that a spin Hall current along the y di¬ 
rection and polarized in the 2 ; direction may exist when 
an external electric field is applied along the x direc¬ 
tion. The symbols b and A are defined as b = h/m 
and A = jdh^, respectively. In the weak anisotropy limit 
[j3 •€. {2mhkF)~^], we keep the leading-order contribu¬ 
tion to the intrinsic SHC and find 


a 


yx 


e 

Stt ’ 


( 21 ) 


which is identical to that of the 2DEG with fc-linear 
Rashba SOC [5S] but is different from 2DEG with the 
/c-cubed Rashba SOC [HI 0^. The vanishment of the 
other components of the spin conductivity tensor is due 
to the symmetry of the Fermi surface. 

Now we consider the impact of disorder on the SHC 
up to the vertex correction. It is more convenient to 
implement the calculation in the eigenvector space. We 
consider the randomly distributed, identical point defects 
that are spin independent, V (r) = aoVo S {r — Ri), 
and the matrix element can be expressed as 




( 22 ) 


where Vq is the strength of defect potential and Ri is the 
position of the defect. The self-energy in the first-order 
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where n = N /is the density of impurities per unit 
area and ((• • • ))av denotes the ensemble averaging over 
the impurity distribution. We have introduced the rela¬ 
tion of disorder-free Green’s function, (fciSi| Go \k 2 S 2 ) = 
<5siS2'^fcifc2 fffcisi■ Thus the disordered Green’s function 
turns out to be 

{{ks\G\k's'))^y = ^ Sss'bk'k = Gks- (24) 

9k8 

For the ladder diagram correction to the velocity opera¬ 
tor, we have the following iterative equation 

(fe) = <..2 (fc) + E E {{Kl'^K'k))AV 

k' ^ 3:^4 

X Gl {k') Gi {k') {k') , (25) 

where is the corrected velocity operator, si, 2 , 3,4 = ±1- 
It is difficult to solve analytically the above self- 
consistent equation due to the anisotropic dispersion. 
However, by considering the weak SOC limit, i.e., 
ImEfcj,. ^ jdkJ’kJp ~ 2 m~' equation can be approxi¬ 
mately solved by keeping the leading order of j3, where 
kp is the Fermi wave vector and Tikp is the self-energy. 

After lengthy but straightforward calculations, we can 
find the corrected velocity operator (its derivation is pre¬ 
sented in the Supplemental Material EH), 


(fc) = (fc) -I- JmEpcry. 


(26) 


Following the similar procedure in Eq. (16 1 , we can cal¬ 


culate the SHC with the vertex correction in the weak 
anisotropy limit and get 




.0-2 1 _ 
yx\Y 
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IGtt 


(27) 


It can be seen that in the weak anisotropy limit, the 
vertex correction reduces the magnitude of SHC by a 
factor of 2. In fact, this unique feature of SHC under the 
influence of disorder originates from the special form of 
Rashba SOC. Our result is qualitatively consistent with 
the fact that the term ak"^ {k x cr) ■ z would result in a 
nonzero SHC even with the vertex correction [43j . On the 
other hand, the vertex correction of disorder can cause 
the intrinsic SHC of 2DEG with standard /c-linear Rashba 
SOC to vanish identically [H], but does not affect the 
one with fc-cubed Rashba SOC |3S]. Hence, the distinct 
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behaviors of SHC can be used to determine the nature of 
the Rashba SOC at LAO/STO interface. 

In summary, we have developed an effective Hamilto¬ 
nian of the 2DEG at the oxide interface LAO/STO and 
found a different anisotropic Rashba SOC. We have found 
that the static spin susceptibilities are anisotropic and 
dependent on the Fermi energy. We have also demon¬ 
strated that this different Rashba SOC possesses entirely 
different behavior for the SHC under disorder. There¬ 
fore, these unconventional spin responses can be used to 
determine the nature of Rashba SOC in experiments. 

We are grateful to Matthew Daniels for a careful read¬ 
ing of the manuscript. This work is supported by AFOSR 
No. FA9550-12-1-0479 and No. FA9550-14-1-0277. 
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SUPPLEMENTARY MATERIAL FOR” SPIN RESPONSES AND EFFECTIVE HAMILTONIAN FOR THE 
TWO DIMENSIONAL ELECTRON GAS AT OXIDE INTERFACE LaAlOa/SrTiOa” 


In this supplementary material, we provide the detailed derivation of the effective Hamiltonians and the corrected 
velocity operator. 


THE DERIVATION OF THE EFFECTIVE HAMILTONIANS 

In this section, we construct an effective Hamiltonian around E point. We first decompose the TB Hamiltonian 
Utb into two parts 


Htb -H{k = 0) + H{k), (28) 

where H [k = 0) and H (k) stand for the Hamiltonian at the exact T point and its deviation, respectively. In a new 
basis {\d^y f), i), \dyz i), \dxy i ),\dxz t ),\dyz t)}, H {k = 0) and H {k) take the following form 


H{k = 0) 


/ —Ae — 

0 

0 

V 0 


-2t 

-2t 
0 0 
0 0 
0 0 


0 

0 

0 

—Ae — 4t 


0 0 \ 
0 0 
0 0 

-2t 

-2t) 


(29) 


and 


H{k) 


/ tk^ 0 0 

0 tkl 0 

0 0 tkl 

0 —2i"/k,i —2i"fkx 

2i-fky 0 0 

y 2i"fkx 0 0 


0 —2ijky —2i^kx \ 

2i"fky 0 0 

2ijkx 0 0 

tk'^ 0 0 

0 tkl 0 

0 0 ) 


(30) 


It can be seen that in H (k) we just keep the leading order term in k for each component with k'^ = kl + kl. At the 
r point, we can get the two-fold degenerate eigenvalues 


Et = -^- 2t 


Em. = 


E,= 


—AE — Qt + ^ /AE + 2t + ^ 


+ 2e 


— Ae — 6t -h ^ ^ Ae 2t ^ 


+ 2^ 


(31) 

(32) 

(33) 


and their eigenstates 
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i’ti 

V’i 2 

' 0 ml 


' 0 m 2 

0 hl 

0 h 2 


C'l ( 6*1 Mxy, t) + i M i) + \dyz,i)) 
Cl { — Ol \dxy, i-) — i \dxzA) + \dyz,t)) 


1 

71 

1 

— \ (ixz-) T) “t” l^yzj 

C2{e2\dxyA) + i\d XZ 1 'I') l^yzt 4')) 

C2 { — O2 \dxy^ i) ~ ^ \dxzjt) + \dyzA)) ■ 


( ^ “ 1 “ 

(i I dxz 1 t) + IdyzA)) 


where the parameters are given by 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


Et + 2t — ^ Ei, + 2t — ^ 

c'l — -T-, ^2 —-T-, Ol ,2 — 


1 






6*1,2 + 2 


Here the subscripts t, m, b refer to top, middle and bottom bands, respectively. 
{'ipti,'4’t2,4’m.i,'4’m2,'ipbi,'fpb2}, the Hamiltonian TItb can be expressed as 


(40) 


In the basis 



( tiifc^ 

iti2k- ti3 

{kl - kl) 

iii4 A)_|_ 

ti3p 

itiQk- 



tiifc^ 

-itiik- 

ti3 {kl - kl) 

—itiek+ 

<15^^ 

Htb = 



tk'^/2 

0 

635 {kl - kl) 

'itsek^. 




tk'^/2 

—iH^k- <35 

{kl - k] 



* 



<55 

itbek- 


1 





hbk'^ 

where the parameters are tn = 

C| (1 + 01] 

) t, ti2 = 

-4C|76»i, ti3 

= —Ci</-\/2, <56 

= -40 


614 — “SCi^iy/v^, 615 — 1 ^ 11^2 (1 + ^ 11 ^ 2 ) t, tiQ — — 2 jCiC 2 {01 + 02 ), 636 — — 2 C 202 JI and tss — C'l (l + 0 |) t. 

We follow the standard quasi-degenerate perturbation theory and seperate the Hamiltonian into the unperturbated 
and perturbated part Htb = 77np + 77per- The corresponding matrix elements of perturbation Hamiltonian i7per read 


, , I I , , _ 6?3 (kl -kl)^ + tl^k^ tl,k^ + tl^k^ 

[Iptl I Tiper mi) - - ^^ -^- ^^ - 

2iti3tu {kl - kl) k++ 2iti5tiek^k. 

( 0 tl| i7per m 2 ) = -^^ 


-‘-'m -‘-'t -*^0 

I , 6 I 3 (fc| — klY + tl^k"^ t§5 {k"l — k^Y + tlgh 

( 0 ml| i7per |^ml) = ^ \ ^ ^ 

Hjt — -C/'W H/h — H/m 


-^riL ^b 

(0ml| i7per |0m2) = 2^ 

, ,, I I ,, 1 _ i?3 {kl -kl)^ + fuk^ tl,k^ + tl^e 

— ^b -^m — ^b 

, I I , 2iti3tiePk- 2it35t36 (fc| - kf) k+ 

(l/^bil i7per |0h2) = - 

— ^b — ^b 

and the others can be obtained 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 




( 0 al| 77per | 0 al) — (' 0 a 2 | TIper |' 0 a 2 ) (^'^) 

(0a2| Hper |0ol)* = (0al | 77per \'lpa2) , (48) 


where a = t,m,b are the band indices. Up to the leading order of the Rashba SOC, we could get the effective 
Hamiltonian for the top pair of bands 



Ht (fc) = 


-iatk+ 2 ^ 


for the middle pair of bands 


Hm (fc) = 


2 ^ ia^ {kl - kl) k+ 


-iam {kl - kl) fc_ 


2 m„ 


and for the bottom pair of bands 


(49) 


(50) 


2^^, i^bk- 
-iabk+ ^ 


Hb (fc) = 

where the effective masses nia and Rashba SOC strength (a = t, m, b) are given by 


(51) 


^ '■M I he 

2mt Ejn — Et Eb — Et 


1 


t 


14 


t 


2m„ 


t 


2 Et — Em Eb — E„ 


36 


t 


16 


+ 


t 


36 


2 m 6 Et-Eb ' Em-Eb 


(52) 

(53) 

(54) 


Cl* — tl2, Otb — ^56 

tistu , ^ 35^36 


am = 2i 


Et — En 


Eb — Ej, 


(55) 

(56) 


THE DERIVATION OF THE CORRECTED VELOCITY OPERATOR 


In this section, we turn to give a detailed calculation of vertex correction. As shown in the main text, we have the 
iterative equation 


(fe) = (fc)+E E {{va'^K':^))Av 

k' S3,S4 

XGS(fc')G^3(fc') ^33.34 (fc'). 


(57) 


where Green’s function = l/{Ep — e^s T with ImE^s = h/2Tks and Si_ 2 , 3,4 = il- Eigenvalue reads 

h'^k'^ 


Sks = -A -H sf3h^k^\ cos26»fc| 

2m 


(58) 


and wave function reads 

ik r / _A /- p-iSk \ 

<(>fc 3 = j, Cfe = cos26<fe/|cos26»fc|. 

Velocity matrix element in the eigenvector space reads 

/ iv%)++ {vl)+-\^f “ cos6»fc + ^^Cfe(cos36»fc + 5cos6»fc) Cfc(sin36>fc - sin6>fc) 

\ivk)-+ (^fc )—) Y ^^.(sin36lfc - sin6>fc) ^ cos0^ - ^^ECfc(cos36»fc + 5cos6»fc) 


(59) 


( 60 ) 
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where {v^)ss' = {4’ks\v^\(t>ks')■ Based on these, one can evaluate the disorder-averaged correlation function 

((htfh'l)) 


= ((htthtJ)),, = fid+- V)). 

nV^ 

° (1-CfcCfc' cos(6»fc 


= ((v^,~ V^~, , , - 

\\ kM k.kllj^y \\ k,k k .k j j jyy \\ k ,k k,k j j 2L'^ 


= - ((n+fn:.-)),,, = - - ey. 


(61) 

(62) 

(63) 


Now we focus on the weak scattering and weak anisotropy limit, i.e., ImSfe^,, <C Ph^kp <C ^nf • sense, we have 

the approximated expression of Fermi wave vector and density of states 

yJ2mEp 2/3m'^Ep 


ks,F 

Ns,p 


m 


h h 

2 


cos 26>fc|, 


3/3m"^ /-——— 
^2-^^V2mEp, 


(64) 

(65) 


where k±^p and N±^p are the Fermi wave vector and density of states for two branches of bands. And the relaxation 
time Tks is given by 


1 


1 




Tfe, 


= IE ((C* Af)),, + T E 


nV^m 
h? ■ 


According to Eq. (60), it is natural to assume that the modified velocity takes the following form 

ivl)++ = (/3AoCfc + Aifc)cos6>fc -l-/3A:^Cfc(^2Cos36»fc -P A3Cos6>fc), 

{Vk)+- = /3Bo(kSm6k + (3k‘^(k{Bism36k + B2sm9k), 

{Vk)-+ = /3CoCkSm0k + Pk'^CkiCi sinS^fc -h C'2sin6>fc), 

(vk) — = {/3DoCk + Dik) cos 9k + fik'^Ck{D 2 cos 36k + cos 9k )■ 


Substitute these matrix elements into Eq. (57), and by using the relations 

kdk 


-^(^k,+^k,- - 


dck- 


iS{Ep - Cfe ) 


2 (f - 3/3n^k\cos29k\) Ep - Ck,- - 2ph^k^\cos29k\ 
f dek,+ iS{Ep - ek,+) 

2(§ + 3/3n^k\ cos29k\) Ep - ek,+ + 2/3;i3p| cos 20 fe| 
m i 


h? 2(3H^kp \ cos 29k \ ’ 
i 


f kdk~p _m 

J 2p '=’+ h'^2f3n^k^p\cos26k\' 

one finally obtains these coefficients of Ai, Bi, Ci and Di, 

A A ^ A A 

Aq = mEp, Ai = —, A 2 = , A 3 = , 

m 2 2 

ifi^ ih^ 

Bo = imEp, Bi = —, B 2 = 

ih'^ ih^ 

Co = -imEp, Ci = - —, C 2 = —, 

h ^h‘2 

Do = -mEp, bDi = -, D 2 = --, Do = - — . 

m2 2 

This indicates a relation between the unmodified and modified velocities (in eigenvector space) 


{v%)++ {vl)+. «)++ 


ivl)-+ m- 


K)-+ «)- 


fdmEpCk 


cos 9k i sin 6k 
—i sin 9k — cos 0 fc 


( 66 ) 


(67) 

( 68 ) 

(69) 

(70) 


(71) 

(72) 


(73) 


(74) 
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When transformed into the spin space, this means 

vl = vl+ PmEpCTy, (75) 

which is nothing but the corrected velocity operator in Eq. (26) in the main text. 



